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1 Introduction and summary 



In this note we address the question of how to understand the spectrum of pure Yang-Mills 
(YM) theory using the formalism of wave-functionals. In particular we review and clarify our 
previous work on 2 + 1 YM theory [1] and discuss some preliminary results regarding the 3 + 1 
dimensional problem. 

The wave functional approach we have recently used in the context of 2 + 1 YM theory pQ 
turns out to fit very nicely the existing lattice simulations [2J. After reviewing this work we 
discuss a particularly simple quasi- Gaussian ansatz for the vacuum wave-functional which can 
be used for a study of the spectrum of 3 + 1 Yang-Mills theory. This ansatz is based on some 
heuristic postulates motivated by existing lattice simulations. Note that even though we cannot 
argue from first principles that we understand the large N limit, we do compare our results to 
the large N lattice data. It turns out that the comparison with finite N results does not alter 
the final results too much. We follow the general strategy towards a wave-functional approach 
to the spectrum of 3 + 1 YM that was outlined in [3]. The lattice results we compare to are 
collected, for example, in 01 El El El • 

The main physics issues to be addressed in 3 + 1 YM theory, as opposed to its 2 + 1 
dimensional counterpart are: 

1) dynamical generation of a scale independent mass: the so called gap equation. We discuss 
this issue by analyzing the origin of the mass gap from the geometry of the configuration space 
in the context of 2 + 1 YM. We then present some very heuristic arguments about the gap 
origin in the more difficult 3 + 1 dimensional problem. 

2) the IR asymptotics of the quasi-Gaussian vacuum wave functional and the issue of the 
evaluation of the string tension of large planar Wilson loops. 

3) finally, there is the question of how to recover the main feature of the UV physics, i.e. 
asymptotic freedom, the hallmark universal feature of QCD perturbation theory (as exemplified 
in the famous 11/3 one-loop (3 function coefficient). This turns out to follow from the proper 
UV asymptotics of the vacuum wave-functional. 

By contrast, one does not have to worry about these issues in the analysis of 2 + 1 YM. 
In that case the theory has a natural dimensionful parameter, the YM coupling. So, in our 
attempt to treat some phenomenological features of the 3 + 1 YM spectrum by using the 
experience with 2 + 1 YM theory, we implicitly assume that the above difficult issues can be 
addressed in a self-consistent manner in the context of simple quasi-Gaussian vacuum wave 
functionals. For example, we simply assume that the mass gap is self-consistently generated in 
the wave functional approach to the 3 + 1 YM theory. Also, the form of the quasi-Gaussian 
ansatz for the vacuum wave functional in the 3 + 1 dimensional context is essentially motivated 
by the 2 + 1 dimensional counterpart. The IR asymptotics of the 3 + 1 YM theory is very 
difficult to analyze from first principles: see, for example, the discussion in [TDj. In the 3 + 1 
problem, the asymptotic freedom issue is dealt with by appealing to the correct UV limit of the 
vacuum wave functional, the requirement of gauge symmetry (so that all interaction vertices 
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are correctly captured) and the standard RG analysis of wave functionals [TT] . 

Our appeal for a self-consistent generation of the mass gap and the IR asymptotics on the 
level of the vacuum wave functional does have some resonance with the philosophy pursued in 
[12J as well as the AdS/QCD approach [13]. The details appear to be very different though. 
The essential physics ingredient in our case is the constituent picture which seems to work so 
nicely in 2 + 1 YM [T] (for previous attempts at the constituent picture in the 3 + 1 YM consult, 
for example, [H] ) 

This note is organized as follows: 

1) We start with a general discussion of how to understand the spectrum of pure Yang Mills 
theory in the formalism of vacuum wave functionals. We concentrate on here 2 + 1 dimensional 
YM theory. The crucial heuristic points in this approach are the issues of a dynamical mass 
gap generation and the role of the kinetic term. 

2) We then present a more phenomenologically minded discussion of the spectrum of 3 + 

1 dimensional Yang Mills theory based on a quasi-Gaussian vacuum wave functional. This 
phenomenologically minded ansatz is on one hand motivated by the results in 2 + 1 dimensions 
and on the other hand by good numerical comparison with the lattice data. We think that 
the results presented in this paper justify the simple quasi-Gaussian vacuum wave-functional 
approach in spite of a somewhat tentative nature of our theoretical analysis. 

2 Wave functionals and the kinetic term 

Here we outline a general procedure regarding the derivation of a vacuum wave functional for 
pure YM theory, which in the case of 2 + 1 YM theory happens to be in excellent agreement 
with the corresponding lattice data concerning the spectrum of the vacuum states of the theory 

ID- 

We contrast the abelian case, which can be exactly solved, and the more difficult non-abelian 
case. We do not utilize any special variables. The approach is, in principle, very general and 
transparent. For concreteness, in what follows, we concentrate on the 2 + 1 dimensional physics. 
In later sections we will discuss the 3 + 1 dimensional YM theory. 

Given the gauge theory Hamiltonian 




where Ei and B are the electric and magnetic components of = d^A u — d u A^ + A„], 
we wish to reason out the vacuum wave-functional using the ordinary connection variables A a . 
In the Hamiltonian gauge A = 




We are looking for ty such that H^ = 0, after subtracting the infinite vacuum energy term. 
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In the abelian case the answer is obviously Gaussian. In the non-abelian case this is not so, 
for reasons of non-abelian gauge symmetry. So, in 2+1 dimensions we consider a quasi-Gaussian 
ansatz, which is compatible with the Gauss law 



^o = exp(-c JB(K(L))B) 



(3) 



where c is the appropriate constant l/(2g 2 m), B is the magnetic field and L = 7—5, where A is 
the covariant Laplacian, and m is the appropriate mass parameter. Note that the parameter m 
can be explicitly evaluated using the approach of "corner variables" [15] of Karabali and Nair 
[j~6] . In that case there is a very precise relation (for large N) 



m = 



g*N 
2tt 



(4) 



It suffices to say here that the parameter m is ultimately related to the value of the string 
tension for large Wilson loops. Also, at least formally 



K(L) = £ c n L n 



(5) 



ra=0 



and thus more generally 



E + J B a (K) B a + . . . 



o- 



(6) 



The (divergent) vacuum energy Eq can be isolated and as expected, the leading divergence in 
the UV is cubic. Next, what we need to do is compute the expression that we have labeled by 
TZ and set it to zero. This will constitute an equation for the kernel K. 

Let us consider the various terms in TZ. First, the potential term B a B a clearly contributes 
a fixed constant to TZ. The remainder of 1Z will come from the action of the kinetic energy 
operator Tkn- Given the form of the vacuum wave functional ^0 ~ e p , it is elementary to 
derive 

5^o 



5P 



5A a (z) 5A a (z 



(7) 



and thus 



from which we find 



5A a (z)A b (w) 



5 2 P 



+ 



5P 



5P 



5A a (z)A b (w) 5A a (z) SA b (w) 



5P 5P 1 
TP + ml °, °, I B a B a 

,w 5A a (z) 5A b (w) m 



(9) 



The second term in brackets is easy to compute. Since we want to solve Schrodinger equation 
to quadratic in B a order only 



P ~ / B a K ( — ) B a + . . . 

,m 2 



(10) 
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By evaluating functional derivatives 5P/5A we see that the second term in brackets in (Q is 
equal to 

A „ 2 ( A 



B c 



—K z — 

m 2 \m 2 



B a (11) 



Therefore, we deduce that this term contributes LK 2 (L) into 1Z. This result also holds in the 
pure Abelian case. 

Let us finally proceed to the TP term in 0. In the Abelian case this term does not depend 
on field variables and gives (infinite) vacuum energy only In the non-Abelian case the kinetic 
term plays a very important role. In 2 + 1 dimensions the simplest way to proceed is to assume 
that the kinetic energy operator T acts homogeneously on any local operator-valued function 
of B. This is natural in the holomorphic corner variable approach of Karabali, Kim and Nair. 
This assumption can be somewhat justified by explicit computations involving the action of T 
on the potential term B 2 . In that particular case the action of T would simply produce 2mP, 
since P is quadratic in B. Now, by extrapolating this statement into [1] 

TO n = (2 + n)mO n + ... (12) 

where O n = BL n B, one gets much more information^. Notice that here ellipsis stand for terms 
of higher order in B, but the same mass dimension as O n , and which mix with O n under the 
action of T. Finally, note that this statement in some sense incorporates the role of topology 
(i.e. compactness of the configuration space) in the non-abelian case, and encapsulates the 
importance of the kinetic term in the dynamical origin of the mass gap, which is intuitively 
very much in line with [TT] . 

Returning to our discussion of the Schrodinger equation we note that ffTSI) implies 



TP = -c J B\^c n (2 + n)L n ^B. 



(13) 



It is convenient to write the factor in braces formally as 

1 d 



L K(L) . (14) 



2LdL 

Assembling all of these results, we then find the following Riccati equation [1] 

L d 



K 



0. (15) 



Note that the abelian answer is given by the last algebraic part of the equation, which does 
not have anything to do with the spectrum of the kinetic term, i.e. LK 2 + 1 = 0. This of course, 
has to be reproduced even in the non-abelian situation at UV. Now, the discrete spectrum of 
T is essential for the emergence of the gap. So the first part of the equation [L 2 K(L)} is 



1 A possible calculational justification for this statement might be obtained by following the strategy proposed 
in the recent work by Fukuma et al [18j . 
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crucial. In the IR limit only the factor K + 1 = survives. The derivative term k*-.[K(L)\ 
interpolates between the UV and the IR limits. The other essential ingredient is that there 
is only one normalizable K that solves the whole differential equation, as discussed in what 
follows. 

Finally, note that the whole discussion is in principle valid for any rank of the gauge group. 
Even though we do not have a very good reason to rely on large N, we will perform the actual 
comparison of analytic results for the spectrum with large N lattice data (which are already 
very close to the 577(3) case.) 



2.1 The Gluonic Kernel 

Although the above Riccati equation is non-linear, it is easily transformed into a linear second- 
order equation of the Bessel type, and one finds a general solution of the form 

= i cM^m+y^m (16) 

y/LCA x (Ay/Z)+Y y {4:VL) 

where C is a constant and J n (Y n ) denote the Bessel functions of the first (second) kind. As 
explained in [1] it is remarkable that the only normalizable wave functional is obtained for 
C — > oo, which is also the only case that has both the correct UV behavior appropriate to 
asymptotic freedom, as well as the correct IR behavior appropriate to confinement and a mass 
gap! This solution is of the form 

K{L) = iM±m (17) 



This remarkable formula is reminiscent of similar results in related contexts [12]; here, it en- 
codes information on the spectrum of the theory. We note that this kernel has the following 
asymptotics (where L ~ — /4m 2 ) 

p^O, K^l; p^oo, K-+2m/p (18) 

consistent with confinement and asymptotic freedom, respectively (Note that the argument of 
Bessel functions is imaginary; so instead of J n we have I n Bessel functions.) 

Now using standard Bessel function identities we may expand 

m i +2u ±^ (19) 

where the 7 2jn are the ordered zeros of ^(u) [I]. The inverse kernel is thus 

K-\L) = v/Z^S = 1 + 8L V - 1 2 . (20) 
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Now, if we regard L ~ dd/m 2 , in terms of momentum p we find 

2^ 1 p 2 + M 2 

Here 

M n = (22) 

M n 's can be interpreted as constituents out of which glueball masses are constructed. It is not 
difficult now to find a Fourier transform of inverse kernel K^ 1 {k). By rewriting (I2TT) as 

i 00 / M 2 \ 

^^^st 1 -^) (23) 

we immediately obtain 

K^ 1 (x -y) = 8®(x - 1/) + 2 E ( * (2) (* - v) ~ ^ ^o(M n |x - y|)J (24) 

where i£o(20 is the modified Bessel function of the third kind. At asymptotically large spatial 
separations \x — y\ — > oo this takes the form 

K-\\x - y\) « £(M n )fe- M «l^l. (25) 

4w27r|x — y| n=i 



As shown in [T] we can treat i? as a free field as well and therefore 

(B(x)B(y))~K-\x-y). (26) 
The constituent masses are controlled by the zeros of J^. 

M„ = -Of- (27) 

The glueball masses then follow by evaluating the correlators of the relevant operators that 
probe the glueball states. For example, the mass of the ++ glueball may be probed by the 
operator TrB 2 . We have pQ 

(Tr(B\Tr(B\}~(K-\\x-y\)) 2 . (28) 



'y 

and thus the ++ mass is pQ 

M 0++ =M 1 + M 1 (29) 

A more complete discussion of other glueball masses can be found in pQ. 

Now it is also not difficult to derive expression for the vacuum expectation of the large 
Wilson loop of area A . At large N 

($(C)) = (TrPexp (i A)) ^exp (-N J d 2 xd 2 y (B(x)B{y))\ (30) 



or 

i(S(C)> = exp J d\d*yK-\x - y)\ (31) 

and from this expression we see that leading 5-function gives area law with a string tension 
that fits beautifully the lattice simulations. As for the rest of terms that appear in (|24p . we 
may notice that 

f A d 2 x^ 2 \x-y)-^K {M n \x-y\)^ - as A - oo (32) 

and therefore these terms will give corrections to the area law behavior which vanish for asymp- 
totically large loops. The string tension then follows to be [16] 



- <J* N F M\ 

Given this expression we can easily convert glueball masses (|29|) into the units of the square 
root of the string tension which makes comparison to lattice data straightforward. 

Finally, a couple of comments. Note that the gluonic wave functional leads to an effective 
partition function of a pure 2d YM theory in the IR limit. Also the purely Gaussian vacuum 
wave-functionals suffice for self-consistent calculations of both the glueball masses and the string 
tension. Similarly, note that the J constituents do not appear as asymptotic states. Thus, it 
might be plausible that the glueball constituents are "seeds" for constituent quarks once the 
fermionic degrees of freedom are included. If so, quark confinement would be "seeded" by the 
confinement mechanism in the pure glue sector. 



2.2 Generalizing the kernel equation 

One of the most intriguing results of [1] is the prediction of the ++ glueball mass which agrees 
extremely well with the lattice simulations. The derivation of that result crucially depends 
on equation (12) which leads to nontrivial kernel equation ( 1T5]) . In the absence of a rigorous 
derivation of that equation from first principles, we want to show that equation (12) is well 
motivated phenomenologically. To this end we may consider a modification of equation (12). 
We will also use the following results based on such a modified kernel equation to motivate our 
phenomenologically minded study of the spectrum of 3 + 1 YM theory. 

Suppose that the action of the kinetic term is as follows^ 

TO n = (2 + bn)mO n + ... (34) 

where b is some free parameter. 

2 In principle we could have taken an even more general expression TO n = (a + bn)mO n + ... . This 
modification, however, does not lead to any new conclusions compared to f|34[) . 
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Then the kernel K(L) satisfies the following equation of the Riccati type 

L d 



After changing variables as 



- K - b 2^l K ( L )}+ LK2 + 1 = ( 35 ) 

K = -b— (36) 
where the prime denotes the derivative with respect to L, we arrive at 

which is readily solved in terms of Bessel functions. The general solution is 

y = Cl L-^ J 2/b ^(A/bVL) + c 2 L 1 -^Y 2/b ^ 1 (A/bVL). (38) 
The physically relevant solution will be 

y = L 1 ^J 2/b _ 1 (4/by/L) (39) 
Using the standard recursive formulae 

J„_l(u) + Jv+l(u) = fj.H 

Ju-i{u) - J v +i(u) = 2J' v (u) 
and analogous formulae for Y u , we arrive at the general solution 



(40) 



K(L) = 1 C J 2/b (4/bVL) + Y 2/b (4/bVZ) 

VLCJ 2/b - 1 (4/bVL)+ Yy^A/bVL)' 

Once again, the only solution which is normalizable in UV and IR is given by 

K( L) = (42) 

We see that the argument (and order) of the relevant Bessel functions depends on b, and we 
recover the standard 2 + 1 result if b = 1. Note that in general if TO n = (a + bn)O n we have 

K{L) = a -±= ^Z^j , (43 ) 
1 ' 2 VI J a/b ^(A/bVL) V } 

For a = 2, the kernel has the following asymptotics (where L ~ — p 2 /Am 2 ) 

p^O, p-*oo, K^2m/p (44) 

consistent with confinement and asymptotic freedom, respectively. 
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Now using standard Bessel function identities (see pLj) or more generally 
where the 7„ jn are the ordered zeros of J u (u). The inverse kernel is thus in general 

k-\l) = yr -v-iWfrff . (46) 

Then the constituent masses in this more general case are given by a simple formula, again 
in parallel with pQ 

U. = (47) 

In the case of 2 + 1 YM the ratio of the glueball masses that follows from this constituent mass 
and the corresponding expression for the square root for the string tension best fits the lattice 
data for b = 1, as discussed in the previous subsection. 



2.3 A recapitulation 

The effect of the kinetic term seems to be of utmost physical importance (at least in the 
context of pure 2 + 1 YM theory.) In the abelian case the kinetic term just contributes to 
the vacuum energy and is not essential for the determination of the kernel. In a topologically 
non-trivial case, i.e. non-abelian theory, the effect of the kinetic term is crucial. Heuristically, 
the compactness of the configuration space seems to imply a discrete "spectrum" for the kinetic 
term, and this drastically changes the nature of the kernel and the resulting mass spectrum 
of the theory. The mass spectrum is in turn very nicely captured by a simple quasi-Gaussian 
vacuum wave functional. Presumably the compact U(l) could be also treated in this manner, 
but it is not clear whether a Gaussian ansatz captures the whole physics. It is possible that 
instead of P ~ / BK(L)B one should consider something non-quadratic yet consistent with 
the compactness of the configuration space, say P ~ / cos(ciBK(L)B), where c\ has the 
appropriate dimensions for the cosine to make sense. 



3 On the origin of the mass gap in pure YM theory 

Here we present a summary of an intuitive argument of Karabali, Kim and Nair concerning the 
mass gap in 2 + 1 YM theory [16]. We then discuss possible generalizations of this argument 
for the case of 3 + 1 YM theory. 

In the Karabali-Nair formalism [TB] one obtains that the configuration space measure is 
given by the exponent of the 2d WZW action. This is the volume element which then also 
defines a normalizable inner product of gauge invariant wave functionals 

(l|2)a+i = / dfi(H)e 2cASwzw{H) ^l(H)^ 2 (H) (48) 
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This configuration space measure is crucial in the argument concerning the dynamical origin of 
the mass gap. (The emphasis here is on the term dynamical; that the gap should be proportional 
to the dimensionful coupling constant is obvious on dimensional grounds.) 

First note that the quadratic part of Swzw{H) is given by, in terms of the magnetic field, 

as 

Swzw ~~r I dx 2 B—B + ... (49) 



Air J ~~ V 2 

where the higher order terms in B are dropped in the exact expression. 

Now, the intuitive argument of [16] runs as follows. Denote by AE and AB the root mean 
square fluctuations of the electric and magnetic fields. By utilizing the canonical commutation 
relations [Ef, AA = —i5ij5 ab , one gets AEAB ~ k, where k denotes the momentum. The 
estimate for the energy is then 

\ s g 2 k 2 AB\ , . 

g being the gauge coupling. Note that this would give, by minimizing over AB 2 , AB 2 min ~ g 2 k 
which then implies e ~ k, i.e. no mass gap. This is indeed the correct answer for the Abelian 
theory, but incorrect for its non- Abelian cousin. 

In the non- Abelian case we need to take care of the non-trivial configuration space measure 
(because the physical meaning of the uncertainty relation hinges on implicitly knowing the 
inner product), which we have in 2 + 1 thanks to [TB]. The equation for the Gaussian part of 
this measure implies that AB 2 ~ 7tk 2 /CA for small values of the momenta, that is in the IR 
limit. Thus, as pointed out in [16], even though the energy e is minimized around AB ~ k, 
the probability, which is dominated by the measure factor, is centered around AB 2 ~ k 2 ir/cA- 
Consequently, the energy scales as e ~ g 2 CA/2n + 0{k 2 ) and indicates the presence of a mass 
gap in the spectrum, proportional to g 2 CA- Again, note the crucial role played by the kinetic 
term in this intuitive argument. 

We wish to emphasize the dynamical nature of this heuristic argument, as opposed to 
the final result which could be expected on dimensional grounds. Obviously the dimensional 
argument fails in 3 + 1, so we discuss some options regarding the generalization of the dynamical 
part of the argument. 

For example, one option would be to consider the natural generalization of the above con- 
figuration space measure for 2 + 1 YM in the 3 + 1 dimensional context 

(1|2) 3+1 = J d^B)e 2 ^ s ^n(B)^ 2 (B) (51) 



S 3+ i~- dx 3 B—B + ... (52) 



where to quadratic order in B, by generalizing the quadratic result from 2 + 1 

m f , o -> 1 

7-l ixB v 

The parameter m should represent an RG invariant mass scale and it should be dynamically 
determined from a gap equation, perhaps along the lines of |10| . (This formula is natural form 
the point of view of the formalism of corner variables discussed in 1191 . A different expression 
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for this measure was found in [20 J. This might indicate the sensitivity of the 3 + 1 dimensional 

problem to a choice of variables and the fact that there are no obvious counterparts to the WZW 
functional in three dimensions.) By repeating the above argument based on the uncertainty 
principle we get that in this case the gap is proportional to m, as expected. 

Another option would be to distinguish the 2+1 and 3+1 dimensional problems by appealing 
to a more dynamical role of the potential term in the 3 + 1 dimensional context. In this case 
the intuition about the infrared dynamics of 3 + 1 YM theory based solely on the experience 
with 2 + 1 YM theory would be most probably incomplete. 

Whatever the precise dynamical mechanism of mass generation is in 3 + 1 YM theory, the 
crucial question is to determine the mass gap dynamically (via a gap equation), and relate it 
to the universal features of perturbation theory (i.e. asymptotic freedom and Aq C d)- We will 
not explicitly address these difficult issue here. In the next section, we will simply postulate 
that the wave-functional "knows" about a stable mass parameter m, which is in turn related 
to the mass gap as well as Aqcd- 

4 A preliminary study of the spectrum of 3 + 1 YM 

In this section, we want to approach the spectrum of 3 + 1 YM theory in the manner used in 
the analysis of 2 + 1 YM, as summarized in the preceding part of the paper. Of course, the two 
theories are radically different, as outlined in the introduction. Our discussion is necessarily 
phenomenologically minded, given the difficult issues concerning the dynamical origin of the 
mass gap and the IR asymptotics of the vacuum wave functional. We will be mainly motivated 
by the existing lattice simulations in 3 + 1 YM and the success of the simple quasi-Gaussian 
ansatz for the vacuum wave functional for 2 + 1 YM theory. 

As in the case of 2 + 1 YM we assume the following ansatz for the vacuum wave functional 



with a non-trivial kernel K and L = ^Aj-, where A is the 3d covariant Laplacian. Note, that in 
contrast to 2 + 1 YM here it is crucial to determine the dynamical scale m (via a gap equation), 
which is not explicitly present in the bare Hamiltonian of 3 + 1 YM. We will model the kernel 
K(L) motivated by the discussion from section 2.2 of a generalized kernel equation for the 2 + 1 
dimensional problem. As pointed out in section 2.2, the UV behavior of such a kernel will be 
just what we need in 3 + 1 YM theory 



In that limit the mass m cancels and we are left with the correct UV result known from 
QCD perturbation theory. By applying the standard background field methods to such a UV 
wave functional (following, for example, the nice work of Zarembo in (TTJ) we can readily 




(53) 




(54) 
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recover the universal features of perturbation theory, such as the famous 11/3 coefficient in 
the one-loop beta function. The corrections to this UV result (indicated by in (I53j) ) 

can be organized in powers of the mass m over the powers of momenta and as such they 
encapsulate non-perturbative effects. Given the usual decoupling between the perturbative and 
non-perturbative effects^, such terms will not spoil the standard results of the QCD perturbation 
theory. In essence, what we propose is to assume that such terms can be summed up to give a 
very particular form for the kernel, motivated by our results regarding the 2 + 1 YM theory as 
well as the existing lattice simulations concerning the spectrum of 3 + 1 YM theory. 

We now proceed with a very particular kernel discussed in the section 2.2 for the case of 
6 = 2/3. To motivate this choice for the parameter b we may naively argue as follows: In 2 + 1 
dimensions the magnetic field B that features in the vacuum wave functional is a scalar. In 
3 + 1 dimensions the magnetic field Bi is a vector, so we naturally have to contract this vector 
index in the vacuum wave functional above. By averaging the scalar product of 3 vectors v l 
in 3d over all directions (a fixed direction being defined by the usual Euler angles 9 and (f>) 
one gets ^v 2 , where 4^ comes from the angular integral / cos 2 9 sin OdOdtp. Per 2tt range of the 
4> angle we also get the factor 2/3. Granting its very naive nature, this might be a heuristic 
way to understand why the 2/3 factor seems to be "natural" in the discussion of the vacuum 
wave functional for 3 + 1 YM given its 2 + 1 dimensional counterpart, but we admit that this 
simplistic argument is too loose and therefore only a phenomenological justification is possible 
at present. 

Thus we take the relevant kernel from section 2.2 by choosing b = 2/3 

K-\L) = v^^S. (55) 



The constituent masses are controlled by the zeros of J; 



3- 



M n = ^p (56) 

Note that this is different from the usual 2 + 1 result M n = 72, ^ m , which involves zeros of J 2 ' 
In order to determine the mass gap (or the mass of the lowest lying glueball) consider the ++ 
states which may be probed by the operator TrB^B^. We have 

(Tr Tr (B«B?) y ) ~ (k~\\x - y\)f . (57) 

The correlation function follows by Wick theorem given our knowledge of the quasi-gaussian 
vacuum wave functional, in complete parallel with [I]. 

Now, following the general procedure outlined in section 2 the mass gap, or the mass of ++ 
is given as a sum of two constituents 

M 0++ =M 1 + M 1 (58) 



3 In (f54|) g should be understood as the running coupling constant 5(A) at some very high scale A >> Aqcd- 
The non-perturbative mass parameter m should be related to 17(A) and A via m ~ Ae -8 '"' l h ° 9 . 
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Table 1: Large N glueball masses in YM$ + \. All masses are in units of the square root of the 
string tension. Comments on the origin of the lattice data are presented in Section 4.1 



State 


Lattice, iV — ► oo 


Theory 


Difference in % 


0++ 


3.307 ±0.053 


3.394 


[2.6] 


0++* 


4.67 ±0.29 


4.294 


[8] 


2++ 


4.80 ±0.14 


5.193 


[8] 


2 ++* 


6.17 ±0.21 


6.058 


[2] 


4 ++ 


7.81 ±0.20 


6.924 


[11] 


4++* 


9.88 ±0.36 


7.779 


[21] 


6++ 


9.34 ±0.55 


8.634 


[7.5] 



and 

M 0++ * = M 1 + M 2 (59) 

etc, as well as 

M 2 ++ = M 2 + M 2 (60) 
etc, as in pQ. Other states follow as in pQ. Note that that relevant zeros of J3 are [21] : 

73,1 = 6.380, 73,2 = 9.761, 73,3 = 13.015, 73,4 = 16.223, 73,5 = 19.409 (61) 

Given the phenomenological nature of our discussion, we note that the zeros of the J3 Bessel 
function do seem to provide the best fit to the available lattice data. 

In order to be able to compare our results for glueball masses with actual lattice data, we 
have to re-express them in terms of some other measurable quantity^ like ^fo or Aqcd- In other 
words, we need to establish the relationship between the string tension and the mass parameter 
m similar to eq.( l33l in (2+l)D. At the moment we do not have any definite proposal on the 
possible generalization of that equation to (3±1)D. To circumvent this difficulty we should 
consider the coefficient in (133]) as an extra fit parameter. Surprisingly enough, the best fit value 
for this parameter turns out to be very close to «^|. For this reason in what follows we simply 
use the 2+1 dimensional result ( 1331) to express glueball masses in terms of the square root of 
the string tension. 

We summarize the actual numerical results in a couple of tables. As noted in the introduc- 
tion, our analysis can be done at any rank N, but we do comparisons with the large N data 
only. The large N data are from [U [5]. There are also many more data available for SU(3) and 
SU (8). The large N extrapolation is done by multiplying these data by \/N 2 — 1/N, as we did 
in 2 + 1. The theoretical results are shown as well as the respective difference from the lattice 
data. As one can plainly see that the zeros of J3 work very well. 

4 However, we would like to point out that eqs. ([58l) - (|60ll can readily be used to obtain definite numerical 
values for ratios of glueball masses. 
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4.1 Comments: lattice vs. theory 



The only genuine large- N lattice estimates are available for three states only 



m ++ 



3.307(53) 



(62) 



m 2 ++ 



4.80(14) 



(63) 



m ++* 



6.07(17) 



(64) 



This is from the most recent study by Lucini, Teper and Wenger [I] and is based on the 
extrapolation of results available for N=2,3,4,6,8. These results are consistent with a previous 
study by Lucini and Teper [5J (based on N =2,3,4,5). The only difference is that error bars are 
smaller in this new simulation (about 1.5 to 3% as opposed to 5-8% in Lucini and Teper.) 

The ++ and 2 ++ mass estimates seem to be reliable. They show a smooth approach to 
the continuum limit for any given value of N as well as a smooth extrapolation from finite N 
values to N = oo. 

The ++ * state is more troublesome. So we should compare critically data from other sources 
as well. Other sources are the PhD thesis of H. Meyer [6] (this includes SU(3) and SU(8) data; 
somewhat surprisingly, Meyer's data differs from Teper et al.) and Morningstar and Peardon [7] 
(SU(3)only; this is considered, at the moment, to be the most precise lattice study for N — 3). 
First, all these authors report that the finite volume corrections are still important for this 
state and so its identification as a glueball is questionable (i.e. this state might be a torelon). 
Second, according to Teper et al. m ++* « 2 m ++, which leads to a possibility that this state 
is actually a two-glueball state with zero relative momentum. Third, the extrapolation of finite 
N values to N = oo is not so good (as can be seen, for example, from the inverted slope of the 
large- N fit in (164")) ) which suggests the possibility that different mass eigenstates were used while 
performing the fit to iV = oo. Finally, there is a large discrepancy between Teper's N = oo 
result m ++* ~ 2 m ++ and the SU(3) result by Morningstar and Peardon m ++* = 1.54 m ++ 
or the SU(8) result by H. Meyer m ++. = 1.42 m ++. 

For these reasons we think that ++ * mass estimate by Teper et al. might be doubted. A 
more reasonable estimate of the ++ * mass can be based on the SU(8) value found in Meyer's 
thesis, which upon rescaling leads to our N = oo estimate of ++ * mass 



This is the value we use in Table 1 and it's compatible with Morningstar and Peardon (i.e. 
approximately 1.5 times heavier than ++ mass). However, we would like to point once again 
that this estimate should be treated with caution. Similarly, the rest of the lattice values in 
Table 1 (for 2++*, 4++, 4++*, 6 ++ states) were obtained from SU{3) and SU{8) data in [6] by 
simple rescaling and therefore should also be treated with caution. 



TUq++* 



4.67 ±0.29 



(65) 
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Note that Meyer's thesis contains 377(3) estimates for the masses of ++ ** and o ++ ***. 
However, these seem to be too heavy to be absolutely trusted (they are in the same ballpark 
as Teper's ++ * estimate) and therefore we choose not to include them in our comparison. 

This concludes our comparison in the PC = ++ sector of the theory and we want to 
proceed now to the PC = — + sector. But before we do that we would like to make a few 
general comments on the qualitative properties of glueball mass spectrum in both 2 + 1 and 
3 + 1 dimensions. It has been known for a while [TH [S] that the low mass glueball spectrum 
can be qualitatively understood in terms of local gluon interpolating operators of minimal 
dimension: higher dimensional operators create higher mass states. In particular, in 2 + 1 
dimensions there is a (mass) dimension-6 operator tr(B 3 ) with quantum numbers and this 
is consistent with the lattice result concerning the existence of a light (with m — « 1.5mo++) 
glueball state in the spectrum of pure Yang-Mills theory in 2 + 1 dimensions. On the 
contrary, in 3 + 1 dimensions modern lattice studies seem to strongly rule out the existence of 
light glueball states [j. Qualitatively this can be explained by the fact that dimension-6 
interpolating operator with quantum numbers simply does not exist in 3 + ID. 

What does exist in 3 + ID however is Oe = fabs(E a x E b ) ■ E c operator with h quantum 
numbers B It should be clear now that in complete analogy with [Tj, i.e. by considering an 
equal-time correlator of two such operators (Oe(x) Oe(v)), we should think of glueball mass 
spectrum in PC = — h sector as being constructed out of three constituent masses. More 
specifically, we obtain 



A comparison of these predictions with lattice data is presented in Table 2. No genuine 
large N lattice data are available for PC = — h states. Therefore, as before, our large N 
"lattice" estimates were obtained by a simple rescaling of the actual SU(3) and SU(8) lattice 
data presented in [6]. 

Finally, we include two comments indicating potential difficulties in our scheme: 

1. As can be seen from Table 2 the seeming large discrepancy for the masses of _+ * and 
2~ + glueballs can be avoided if we swap theoretical values for these states. After all, there is no 
particular reason why Mi + Mi + M 2 combination should be associated with CT + * and not with 

5 The only definite estimate of 0~~ mass, m -- = (2.44±0.25)to ++ , can be found in [5] and is indeed quite 
heavy. 

6 In fact, in 3 + ID it is possible to construct also a dimension-4 interpolating operator with h quantum 
numbers, which is E a ■ B a . Similarly, one can construct ++ operators (E a ) 2 — (B a ) 2 and f a bc(B a x £> ) ■ B c of 
dimension four and six respectively. However as was argued previously in [14] , for a broad class of constituent 
glue models the lightest h glueball couples to dimension-6 operator only. This explains our choice of Oe = 
fabs(E a x E b ) ■ E c as a probe operator for this state. 



M -+ = 

M 2 -+ 
M 2 -+* 



Mi + Mi + Mi = 5.091 y/Z 
-- Mi + Mi + M 2 = 5.99 y/a 
-- Mi + M 2 + M 2 = 6.89 yfa 
.Ui • M-> • .U :! 7.76 y/a 



(66) 
(67) 
(68) 
(69) 
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Table 2: Large N glueball masses in YM 3+ i. All masses are in units of the square root of the 
string tension. Comments on the origin of the lattice data are presented in Section 4.1. 



State 


Lattice, iV — > oo 


Theory 


Difference in % 


o-+ 


4.81 ±0.13 


5.091 


[6] 


-+* 


7.22 ±0.33 


5.99 


[22] 


2-+ 


5.95 ±0.10 


6.89 


[15] 


2 -+* 


7.46 ± 0.29 


7.76 


[4] 



2~ + . Same for Mi ± M 2 ± M 2 combination. On the other side this is similar to the situation 
involving ++ * and 2 ++ glueballs: here the lattice data favor m ++* > m 2 ++. (However, the 
final status of the ++ * mass on the lattice is far from clear.) Nevertheless, our model predicts 
otherwise. 

2. The 0~ + mass agrees nicely with the lattice result. However, there is a claim in the 
lattice literature that h should be heavier than 2 ++ . This can not be seen from large-N 
extrapolated data in Tables 1 and 2 (they seem to be simply degenerate) but Morningstar and 
Peardon say that "the pseudoscalar is clearly resolved at 7a level to be heavier than the tensor" . 
Our current theoretical model (based on zeros of J3) predicts just the opposite (compare Tables 
1 and 2). It's easy to see however that M -+ > M 2 ++ requirement translates into 7^2 < § lv,\ 
requirement for Bessel zeros. This is possible if v > 3.3. Therefore, from this point of view 
v = 3 is not completely satisfactory! 

In spite of these important caveats we take the above phenomenological fit to the lattice 
data as a positive indication that the general form of the wave functionals we have discussed in 
this note does capture the essential non-perturbative physics. Obviously, we are only scratching 
the surface in our understanding of what is truly going on here. Much more detailed work is 
needed to elucidate the real nature of the wave functionals discussed in this paper. 
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